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We investigate the critical behavior of continuous (second-order) phase transitions in the context of 
(2+l)-dimensional Ginzburg-Landau models with a double-well effective potential. In particular, we 
show that the recently-proposed configurational entropy (CE)-a measure of the spatial complexity 
of the order parameter based on its Fourier-mode decomposition-can be used to identify the critical 
point. We compute the CE for different temperatures and show that large spatial fluctuations near 
the critical point (T c )-characterized by a divergent correlation length-lead to a sharp decrease in 
the associated configurational entropy. We further show that the CE density has a marked scaling 
behavior near criticality, with the same power of Kolmogorov turbulence. We reproduce the behavior 
of the CE at criticality with a percolating many-bubble model. 
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INTRODUCTION 


From materials science pQ to the early universe |2], 
phase transitions offer a striking illustration of how 
changing conditions can affect the physical properties of 
matter [3j. In very broad terms, and for the simplest 
systems described by a single order parameter, it is cus¬ 
tomary to classify phase transitions as being either dis¬ 
continuous or continuous, or as first or second order, re¬ 
spectively. First-order phase transitions can be described 
by an effective free-energy functional (or an effective po¬ 
tential in the language of field theory) where an energy 
barrier separates two or more phases available to the sys¬ 
tem. The system may transition from a higher to a lower 
free-energy state (or from a higher to a lower vacuum 
state) either by a thermal fluctuation of sufficient size (a 
critical bubble) or, for low-temperatures, by a quantum 
fluctuation. Generally speaking, this description of ther¬ 
mal bubble nucleation is valid as long as fF^^/ksT 1 , 
where fF[(j)b\ is the 3d Euclidean action of the spherically- 
symmetric critical bubble or bounce 0&(r), ks is Boltz¬ 
mann’s constant, and T is the environmental tempera¬ 
ture. For quantum tunneling, one uses instead S^[4>b]/h, 
where £4 [00 is the 0(4)-symmetric Euclidean action of 
the 4d bounce. 

For second-order transitions the order parameter varies 
continuously as an external parameter such as the tem¬ 
perature is changed [3 . A well-known example is that 
of an Ising ferromagnet, where the net magnetization 
of a sample is zero above a critical temperature T c , the 
Curie point, and non-zero below it. Below T c the tran¬ 
sition unfolds via spinodal decomposition, whereby long- 
wavelength fluctuations become exponentially-unstable 
to growth. This growth is characterized by the appear¬ 
ance of domains with the same net magnetization which 
compete for dominance with their neighbors. In the 
continuum limit, systems in the Ising universality class 
can be modeled by a Ginzburg-Landau (GL) free-energy 
functional with an order parameter 0(x) [4]. In the ab¬ 


sence of an external source (or a magnetic field), the GL 
free-energy functional is simply 






( 1 ) 


where t = (T — T c ), and 7, a, b are positive constants. For 
T > T c the system has a single free-energy minimum at 
0 = 0, while for T <T C there are two degenerate minima 
at 0o = =L(—at/ 6 ) 1 / 2 . This mean-field theory description 
works well away from the critical point. In the neighbor¬ 
hood of T c one uses perturbation theory and the renor¬ 
malization group to account for the divergent behavior 
of the system. This behavior can be seen through the 
two-point correlation function G( r): away from T c G( r) 
behaves as ~ exp[—r/£(T)], where £(T) is the correla¬ 
tion length, a measure of the spatial extent of correlated 
fluctuations of the order parameter. In mean-field the¬ 
ory, £(T) ~ \T — T c \~ v , where v = 1/2, independently of 
spatial dimensionality. In the neighborhood of T c , where 
the mean-field description breaks down, the behavior of 
spatial fluctuations is corrected using the renormaliza¬ 
tion group. Within the e-expansion, one obtains, in 3d, 
vw 1/2 + e/12 # 7e 2 /162 ~ 0.63 0. 

For continuous transitions in GL-systems, the focus 
of the present letter, the critical point is characterized 
by having fluctuations on all spatial scales. This means 
that while away from T c large spatial fluctuations are 
suppressed, near T c they dominate over smaller ones. 
In this letter, we will explore this fact to obtain a new 
measure of the critical point based on the information 
content carried by fluctuations at different momentum 
scales. For this purpose we will use a measure of spatial 
complexity known as configurational entropy (CE), re¬ 
cently proposed by Gleiser and Stamatopoulos [5]. The 
CE has been used to characterize the information con¬ 
tent [5 and stability of solitonic solutions of field theo¬ 
ries |3 E], to obtain the Chandrasekhar limit of compact 
Newtonian stars [ 6 |, the emergence of nonperturbative 
configurations in the context of spontaneous symmetry 
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breaking [7] and in post-inflationary reheating [ 8 ] in the 
context of a top-hill model of inflation (J)]. (Note that our 
usage of the name “configurational entropy” differs from 
other instances in the literature, as for example in protein 
folding mao Here we show that in the context of con¬ 
tinuous phase transitions the CE provides a very precise 
signature and a marked scaling behavior at criticality. 

MODEL AND NUMERICAL IMPLEMENTATION 

We consider a ( 2 + 1 )-dimensional GL model where the 
system is in contact with an ideal thermal bath at tem¬ 
perature T. The role of the bath is to drive the sys¬ 
tem into thermal equilibrium. This can be simulated 
as a temperature-independent GL functional (so, setting 
t = — 1 in Eq. [lj with a stochastic Langevin equation 

A + r]X — V 2 A — A + A 3 + £ = 0, ( 2 ) 

where we have introduced dimensionless variables = 
and = y/f-A, and we took 7 = 1. In Eq. 
[ 2 J 77 is the viscosity and £(x, £) is the stochastic driv¬ 
ing noise of zero mean, (£) = 0. The two are related 
by the fluctuation-dissipation relation, (£(x, t)£(x', t')) = 
2 r]0S(x — x.')S(t — t'). 0 = T/a is the dimensionless tem¬ 
perature and we take ks = h = 1 . 

We use a staggered leapfrog method and periodic 
boundary conditions to implement the simulation in a 
square lattice of size L and spacing h. We used L = 100 
and h = 0.25. Simulations with larger values of L pro¬ 
duce essentially similar results, apart from typical finite- 
size scaling effects [4]. Different values of h can be renor¬ 
malized with the addition of proper counter terms, as 
has been discussed in Refs. E3 and m- Since all we 



FIG. 1: (Color online.) The order parameter (A) [top (blue) 
line] and the volume fraction (pv) occupied by the A > 0 
phase [bottom (green) line] vs. temperature. The shaded 
regions correspond to la deviations from the mean. Within 
the accuracy of our simulation, the critical temperature is 
0 C ~ 0.43 ± .01, marked by the vertical band. 


need here is to simulate a phase transition with enough 
accuracy, we leave such technical issues of lattice im¬ 
plementation of effective field theories aside. We follow 
Ref. [12 for the implementation of the stochastic dy¬ 
namics, so that the noise is drawn from a unit Gaussian 
scaled by a temperature dependent standard deviation, 
£ = ^2r]6/Ath 2 . To satisfy the Courant condition for 
stable evolution we used a time-spacing of At = h/ 4. 
The discrete Laplacian is implemented via a maximally 
rotationally invariant convolution kernel [15] with error 
of 0(h 2 ). 

We start the field A at the minimum at A = — 1 and 
the bath at low temperature, 0 = 0.01. We wait until the 
field equilibrates, checked using the equipartition theo¬ 
rem: in equilibrium, the average kinetic energy per de¬ 
gree of freedom of the lattice field, (A z 2 /2), is (X 2 -) = 0. 
Once the field is thermalized, which typically takes about 
1,000 time steps, we use ergodicity to take 200 readings 
separated by 50 time steps each to construct an ensemble 
average. We then increase the temperature in increments 
of 0.01 and repeat the entire procedure until we cover the 
interval 6 G [0.01,0.83]. 

We then obtain the ensemble-averaged (A) vs. 0. The 
coupling to the bath induces temperature-dependent fluc¬ 
tuations which, away from the critical point, can be de¬ 
scribed by an effective temperature-dependent potential, 
as in the Hartree approximation m- The critical point 
occurs as (A) —» 0 , when the Z 2 symmetry is restored. 

In Fig. [l]we plot the results. The top (blue) line is (A), 
while the bottom (green) line is the ensemble-averaged 
fraction of the volume occupied by A > 0 (pv)- Symme¬ 
try restoration corresponds to this fraction approaching 
0.5. Shadowed regions correspond to lcr deviation from 
the mean. Within the accuracy of our simulation, the 
critical temperature is 0 C ~ 0.43 ± .01. In the top row 
of Fig. [2] we show the field at different temperatures, in¬ 
cluding at ~ T c , where large-size fluctuating domains are 
apparent, indicative of the divergent correlation length 
at 6 C . 

CONFIGURATIONAL ENTROPY OF THE 
CRITICAL POINT 

Consider the set of square-integrable bounded periodic 
functions with period L in d spatial dimensions, /(x) G 
L 2 (M d ), and their Fourier series decomposition, /(x) = 
£ kn F(k„)e* k -- with k„ = 2n(ni/L,...,n d /L), and 
rii integers. Now define the modal fraction /k n = 
|F(k n )| 2 / |F(k n )| 2 . (For details and the extension to 

nonperiodic functions see [5].) The configurational en¬ 
tropy for the function /(x), Sc[f], is defined as 

S c[f] = - F /kn M/k„]- (3) 

The quantity cr(k n ) = —/k n l n [/k n ] gives the relative en- 
tropic contribution of mode k n . In the spirit of Shan- 
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FIG. 2: (Color online.) Top Row: Snapshots of the equilibrium field X(x,y) for different temperatures. The bar on the right 
denotes the field magnitude. Bottom Row: The corresponding mode distribution of the CE density. Far from the critical 
temperature, the CE density is scale-invariant for low k , while close to the critical temperature flow into IR produces scaling 
suggestive of turbulent behavior, with power-law k~ 5 ^ 3 as in Kolmogorov turbulence m- 


non’s information entropy HU, s c [/] gives an informa¬ 
tional measure of the relative weights of different k- modes 
composing the configuration: it is maximized when all 
N modes carry the same weight, the mode equipartition 
limit, /k n = l/N for any k n , with Sc[f] = In TV. If only 
a single mode is present, Sc[f] = 0. For the lattice used 
here, with N = 400 2 points, the maximum entropy is 

£max = 1L98> 



Temperature 

FIG. 3: (Color online.) Configurational entropy vs. 0 for the 
ensemble-averaged field ( X ). The minimum at 0 C is apparent. 
The lighter shade of grey corresponds to a lcr deviation from 
the mean. 

Consider what a field profile looks like for the above 
examples. Plane waves in momentum space have equally 
distributed modal fractions, and their position space rep¬ 
resentations are highly localized. Conversely, singular 
modes in momentum space have plane wave representa¬ 
tions in position space which are maximally delocalized. 
Localized distributions in position space maximize CE 


(many momentum modes contribute), while delocalized 
distributions minimize it. Sc[f] is, in a sense, an entropy 
of shape, an informational measure of the complexity of 
a given spatial profile in terms of its momentum modes. 
The lower Sc[f], the less information (in terms of con¬ 
tributing momentum modes) is needed to characterize 
the shape. In the context of phase transitions, we should 
expect Sc[f] to vary at different temperatures, as dif¬ 
ferent modes become active. In particular, given that 
near T c the average field distribution is dominated by 
a few long-wavelength modes, we should expect a sharp 
decrease in the CE. 

We use Eq. [3] to compute the configuration entropy of 
(j)(x,y). The discrete Fourier transform is 

&(k x ,ky) = Af 52 4>(n x h,n y h)e~ ih{n * K+n v k v\ 

n x = 0 n y = 0 

(4) 

where k x ^ = 2'Km x ( y )/L, and m x ( y ) are integers. We 
display the CE density as a function of mode magnitude 
\k\ in the bottom row of Fig. [2] For temperatures far 
from the critical value, we find a scale-invariant CE den¬ 
sity for long wavelength modes, and a tapering distribu¬ 
tion for short wavelengths. As the critical temperature 
is approached, the field configurations are dominated by 
large wavelength modes, and the CE density flows toward 
lower values of |fc|, disrupting the scale invariance of the 
spectrum. This is expected from the diverging correla¬ 
tion length that characterizes the critical point. We find 
a scaling behavior near T c with slope & -5 / 3 , the same 
as Kolmogorov turbulence in fluids m ■ Following the 
analogy, we call this critical scaling information- entropic 
turbulence , here with modes flowing from the UV into 
the IR. 

The flow of CE-density into IR modes illustrates how 
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lower -k modes occupy a progressively larger volume of 
momentum space as the critical point is approached. 
Hence the sharp decline in CE near T c . In fact, a simple 
mean-field estimate predicts that the CE will approach 
zero at T c in the infinite-volume limit: using that in 
mean-field theory £(T) ~ \T— T c | -1 / 2 and considering the 
dominant fluctuations far away from T c as being Gaus- 
sians with radius £(T), we can use the continuum limit 
to compute the CE of a single Gaussian in d = 2 as |5j: 
Sc[T] = 27r/^ 2 = 47r(l — T/T c ). Although even far away 
from T c we shouldn’t expect this simple approximation 
to match the behavior of Fig. [3J the general trend is 
for CE to vanish at T c in the infinite volume limit. In 
a finite lattice of length L, there is going to be a mini¬ 
mum value for CE given by the largest average fluctua¬ 
tion within the lattice. Indeed, in Fig. [ 3 ] we see that the 
critical point is characterized by a minimum of the CE 
at Sc{O c ) ~ 0.80 =b 0.30. We found that as a function of 
\0 — 6 c \ the CE drops super-exponentially as criticality is 
approached from below, whereas from above we obtain 
the approximate scaling behavior Sc[4>\(0) oc \0 — 6 c \~^. 

For a more realistic estimate of the minimum value of 
CE for a finite lattice near criticality, we model a large 
fluctuation as a domain of radius R and a kink-like func¬ 
tional profile 


4> 0 , v) 


1 

2 


1 + tanh 


~ a(<p)r ~y 


( 5 ) 


where d measures the thickness of the domain wall in 
units of lattice length L and a is a random perturbation 
defined as 


10 

a((f) = l + V'— cos(n^ + /?„), (6) 

“ n 

n =3 

with a n and /3 n being uniformly-distributed random 
numbers defined in the intervals (0,1) and (0,27r), re¬ 
spectively. We took d = 0.01T so that the domain 
wall thickness matches the zero-temperature correlation 
length. Eqs. [5] and [6] define a fluctuation interpolating 
between 0 = 0 and near 0 = 1. In Fig. [4] we show 
the results for an ensemble of such fluctuations occupy¬ 
ing different volume fractions. Note how the symmetric 
fluctuation (dashed line, obtained by setting a(ip) = 1 in 
Eq. [5| has lower CE than the asymmetric fluctuations, as 
one should expect given that CE measures spatial com¬ 
plexity. Fluctuations with R > L/ 4 (volume fraction 
> 0.2) begin to have boundary issues due to the tail 
of the configuration. If we thus take a fluctuation with 
R/L = 0.25 to represent large fluctuations near T c , from 
Fig. [4] its ensemble-averaged CE is CE ~ 3.72 zb 0.13. Of 
course, a single bubble doesn’t match the complexity of 
the percolating behavior at T c . To simulate the perco¬ 
lating transition, we placed different numbers of bubbles 
randomly on the lattice with initial radius R/L = 0.001. 
We then let their radius R grow to R/L = 0.1, measuring 


the volume fraction (py) they occupy, while computing 
the CE as their radii increase. In Fig. [5] we show the 
CE as a function of py. From bottom up, the number 
of bubbles is 5, 10, 25, 50, 100, and 150. The dashed 
line corresponds to a symmetric bubble. As the number 
of bubbles increases, the CE at py = 0.5 approaches the 
numerical value at T c (see Fig. [3|. The results show 
a simple behavior, Sc(N,py) ~ B(N) log(py/L 2 ), with 
B(N) having a weak N dependence. 

We have presented a new diagnostic tool to study 
critical phenomena based on the configurational entropy 
(CE). We have shown how criticality implies in a sharp 
minimum of CE, and identified a transition from scale- 
free to scaling behavior at criticality, with the same 
power law as Kolmogorov’s turbulence. We introduced 
a percolating-bubble model to describe our numerical re¬ 
sults. We are currently exploring the notion of informa¬ 
tional turbulence and its relation to more complex per¬ 
colation models with varying bubble sizes and hope to 
report on our results soon. 



FIG. 4: (Color online.) Configurational entropy vs. the log 
of volume fraction occupied by a single fluctuation with tanh 
profile defined in Eqs. [5] and [6] The shadowed area corre¬ 
sponds to la deviation from mean. The dashed line corre¬ 
sponds to a symmetric bubble. 
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